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Abstract 

We derive the microcanonical partition function of the ideal relativistic quantum 
gas of spinless bosons in a quantum field framework as an expansion over fixed mul- 
tiplicities. Our calculation generalizes well known expressions in literature in that 
it does not introduce any large volume approximation and it is valid at any volume. 
We discuss the issues concerned with the definition of the microcanonical ensemble 
for a free quantum field at volumes comparable with the Compton wavelength and 
provide a consistent prescription of calculating the microcanonical partition func- 
tion, which is finite at finite volume and yielding the correct thermodynamic limit. 
Besides an immaterial overall factor, the obtained expression turns out to be the 
same as in the non-relativistic multi-particle approach. This work is introductory to 
derive the most general expression of the microcanonical partition function fixing 
the maximal set of observables of the Poincare group. 
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1 Introduction 

The microcanonical ensemble of the relativistic e; cLS IS cL subject which has 
not received much attention in the past. The reason of scarce interest in this 
problem is the peculiarity of physical applications, which have been essentially 
confined within statistical approaches to hadron production and the bag model 
[1]; these are indeed the only cases where involved volumes and particle num- 
bers are so small that microcanonical corrections to average quantities become 
relevant. Otherwise, the involved energies or volumes are so large that canoni- 
cal and grand-canonical ensembles are appropriate for most practical purposes 
(e.g. in relativistic heavy ion collisions). Recently [2], it has been pointed out 
that the equivalence, in the thermodynamic limit, between grand-canonical, 
canonical and microcanonical ensembles does not apply to fluctuations, more 



in general to moments of multiplicity distributions of order > 1. Indeed, effects 
of the difference between statistical ensembles might be unveiled in studying 
multiplicity distributions in relativistic nuclear collisions, [3]. In view of these 
phenomenological applications, it would be then desirable to have an in-depth 
analysis of the microcanonical ensemble of a relativistic quantum gas. 

The main difficulty in tackling this problem stems from the need of imposing a 
finite volume. This is necessary to have a correct thermodynamic limit because, 
being the energy E finite by construction, also V must be finite if the limit 
with E/V fixed is to be taken. Strange as it may seem, a full and rigorous 
treatment of the relativistic microcanonical ensemble of an ideal gas at finite 
volume is still missing. In all previous works on the subject, at some point, the 
large volume approximation is introduced; this is tacitly done, for instance, 
considering the single-particle level density as continuous, namely replacing 
sums over discrete quantum states with momentum space integrations [4]: 



In a previous work [5], we have derived an expression of the microcanoni- 
cal partition function of an ideal relativistic quantum gas with explicit finite- 
volume corrections (see eqs. (23), (24)). However, that expression was obtained 
in an essentially multi-particle first-quantization framework, which, as pointed 
out in ref. [6] should be expected to become inadequate at very low volumes, 
comparable with the Compton wavelength of particles. In this regime, under- 
lying quantum field effects should become important and pair creation due to 
localization an unavoidable effect [7]. Indeed, there are several studies of the 
microcanonical ensemble of a free quantum field [8] , but all of them, again, at 
some point invoke a large- volume approximation. In the limit of large volumes 
one obtains the same expressions of the microcanonical partition function and, 
consequently, of statistical averages as in the first-quantization multi-particle 
approach followed in ref. [5]. 

The aim of this work is to derive a general expression of the microcanonical 
partition function in a full relativistic quantum field framework, valid for any 
finite volume, generalizing the results obtained in [5]. We will do this for the 
simplest case of an ideal gas of spinless bosons and postpone the treatment 
of particles with spin to a forthcoming publication [9]. We will show that 
the expression of the microcanonical partition function obtained in ref. [5] 
in a non-relativistic multi-particle approach holds provided that a consistent 
prescription of subtracting terms arising from field degrees of freedom outside 
the considered volume is introduced. 

The paper is organized as follows: in Sects. 2 and 3 we will argue about general 
features of the microcanonical ensemble for a relativistic system and discuss 
several issues concerning a proper definition of the microcanonical partition 
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function. In Sect. 4 we will cope with the further issues related to the definition 
of a microcanonical ensemble for a quantum field at finite volume. Sects. 5 
and 6 include the main body of this work, where the microcanonical partition 
function is worked out in a quantum field theoretical framework. In Sect. 7 we 
will summarize and discuss the results. 



2 On the definition of the microcanonical partition function 

It is well known that the fundamental tool to calculate statistical averages in 
any ensemble is the partition function. For the microcanonical ensemble one 
has to calculate the microcanonical partition function (MPF) which is usually 
defined as the number of states with a definite value E of total energy: 

ft = ]T 5(E - £ state ) • (2) 
states 

For a quantum system, the MPF is the trace of the operator 5(E — H): 

ft = tid(E - H) (3) 

with proper normalization of the basis states. For instance, for one non- 
relativistic free particle, one has to calculate the trace summing over plane 
waves normalized so as to (p|p') = <5 3 (p — p') : 

n = u ( e -h)^ S (e-Q (p|p) = ^ / A/dV (e - Q 

(4) 

Thereby, one recovers the well known classical expression implying that the 
MPF is the number of phase space cells with size h 3 and given energy E. This 
number is infinite as the volume is unbounded and it is thus impossible to 
calculate a meaningful thermodynamic limit at finite energy density. Hence, 
one usually considers a system confined within a finite region by modifying the 
hamiltonian H with the addition of infinite potential walls, i.e. setting H' = 
H + V where H is the actual internal hamiltonian and V an external potential 
implementing infinite walls. Classically, this leads to a finite ft, namely: 



Throughout this work quantum operators will be distinguished from ordinary 
numbers by a symbol " " " . 
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where V is the volume of the region encompassed by the potential walls 
Also the corresponding quantum problem can be easily solved and one has: 

II = tr(B - A*) S (6) 

where the sum runs over all wave vectors k which, for a rectangular box 
with side Lj and periodic boundary conditions, are labelled by three integers 
(ni,?^, ^3) such that k; = n^jLi. The difficulty of the quantum expression 
(6) with respect to the classical one (5) is that, for a given energy E, a set of 
integers fulfilling the constraint imposed by the Dirac's 5: 

+ + = F 

in general does not exist. Therefore, the MPF vanishes except for a discrete set 
of total energies, for which it is divergent. One has a finite result only for the 
integral number of states J E dEQ(E), that is the number of states with an 
energy less than a given E', but this is clearly a stepwise and non-differentiable 
function of E'. 

This holds for an ideal gas of any finite number of particles: strictly speaking, 
the MPF cannot be defined at finite energy and volume as a continous func- 
tion. Only in the thermodynamic limit E —>■ 00 and V —>■ 00 an expression like 
(6) becomes meaningful because it is then possible to replace the sum over 
discrete levels with a phase space integration: 



cells 



Therefore, for a truly finite quantum system, one needs a better definition 
of microcanonical partition function. A definition which does not suffer from 
previous drawbacks is the following: 

n = tx v 6(E-H) = J2(hv\S{E-H)\h v ) . (8) 

hy 

where H is the internal hamiltonian, without external confining potential, and 
the |foy)'s form a complete set of normalized localized states, i.e. a complete 
set of states for the wavefunctions vanishing out of the region V. It should 
be stressed that these states are not a basis of the full Hilbert space because 
wavefunctions which do not vanish out of V cannot be expanded onto this 
basis; thence the notation try instead of tr meaning that the trace in eq. (8) 
is not a proper one. The difference between (8) and a definition like (3) is 
that \hy) are not eigenstates of the hamiltonian H and the right hand side of 



We will use the same symbol V to denote both the finite region and its volume. 
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(8) does not reduce to a discrete sum of 5's. In fact, this is crucial to have a 
continuous function of E, unlike (6). 



As an example, let us work out the definition (8) for the single free particle 
confined in a rectangular box by infinite potential walls and compare it to (4). 
A complete set of states for this problem is: 

-^exp(ik-x) if x G V k = nn x /L x 7ii + n y / L y jirn z / L z k 
if x i V 

(9) 

Therefore, the MPF definition (8) implies: 



n = £<k|5(£ - £)|k) = E/ d 3 P |(k| P )| 2 5 \e - |y (io) 

where we have inserted a resolution of the identity by using a complete set of 
states for the full Hilbert space. The sum in (10) can be calculated and yields: 



/ d 3 x exp[i(k — p) ■ x] 

J V 

= E 77 I d3x / d3x ' ex Pt lk • ( x - x ')] exp[-ip • (x - x')] 

= (2^X d3x X dV * ,(x - x ' )|exp| - ip ' (x 

where the completeness relation in V: 

E^exp[ik.(x-x')]=% 

k 

has been used. Thus, by using (11), the eq. (10) turns into: 

n = -^- /d 3 P <5(£-— 1 (13) 

(27r) 3 i 1 \ 2mj K ' 

that is the same expression (5) as in the classical case. 

The MPF (13) is now manifestly a continuous function of E and, remarkably, 
its thermodynamic limit V — > oo is the same as the thermodynamic limit of 
the "pure" quantum expression (6) (because of (7)). Since the only strict re- 
quirement for a well defined MPF is to reproduce the correct thermodynamic 
limit, for a gas one can choose a definition like (8) instead of Q = tr8(E — H') 
in (6). We emphasize again that in the passage from (6) to (8) the hamilto- 
nian embodying an external confining potential is replaced with the internal 



E 



(k| P )| 2 = E 



V(2tcY 



v 



x = 



-x) 



(12) 
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hamiltonian while, at the same time, the localized eigenstates of the former 
hamiltonian are used to calculate the trace. 



3 The microcanonical partition function of a relativistic system 

In special relativity, the microcanonical ensemble must include momentum 
conservation beside energy's to fulfill Lorentz invariance. This means that the 
MPF definition (8) should be generalized to [5]: 

Q = ^2(h v \S\P-P)\h v ) , (14) 

hy 

P being the four-momentum of the system and P the four-momentum op- 
erator. The MPF now being a number of states per four-momentum cell, it 
is a Lorentz-invariant quantity. The calculation of the MPF is easiest in the 
rest-frame of the system, where P = (M, 0) and the four-volume Vu, u being 
the four- velocity and V the proper volume of the system, reduces to (V,0), 
according to the usual formulation of statistical relativistic thermodynamics 
[10,11]. 

The eq. (14) can be further generalized by enforcing the conservation of not 
only energy-momentum but of the maximal set of conserved quantities, i.e. 
a maximal set of commuting observables built with the generators of the 
Poincare group. To achieve this, one has to replace 5 A (P — P) in (14) with 
a generic projector Pj over an irreducible state of the representation of the 
Poincare group [5,6], i.e: 

^ = EM p *IM . (is) 

hy 

This ensemble is still generally defined as microcanonical ensemble and (15) 
microcanonical partition function. 

In this work, we will confine to a microcanonical ensemble where only energy- 
momentum is fixed, i.e. our projector P, in eq. (15) will be: 

P l = 5\P-P) (16) 

and to an ideal quantum gas, i.e. with P being the free four-momentum op- 
erator. In fact, it should be stressed that 5 4 (P — P) is not a proper projector, 
because P 2 = aP where a is a divergent constant. This is owing to the fact 
that normalized projectors onto irreducible representations cannot be defined 
for non-compact groups, such as the space-time translation group T(4). Nev- 
ertheless, we will maintain this naming even for non-idempotent operators, 
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relaxing mathematical rigor, because it will be favorable to adopt the projec- 
tion formalism. 

It is worth pointing out that the definition (8) involving only the internal (free) 
hamiltonian, is much more fit than (6) for a relativistic generalization. Besides 
the advantage of restoring continuity in E, discussed in previous section, this 
formulation can be easily extended to the full set of conservation laws with- 
out major conceptual difficulty Conversely, had one tried to generalize (6), 
one should have defined a finite region and afterwards sought the observables 
commuting with the hamiltonian supplemented with an external confining po- 
tential. This would have not been an easy task, and, moreover, a maximal set 
of observables commuting with the modified H would not, in general, define a 
Poincare algebra. This is a well known problem in the static bag model where 
the translational invariance is manifestly broken and momentum is thus not 
conserved. On the other hand, in the definition (15), we deal with the original 
Poincare algebra of unmodified (free) operators and enforce the localization 
through the projector onto confined states. 

The eq. (15) can be recast as a full trace by inserting a complete set of states 
\f) into (15): 

n=EWEI/X/l p *IM 

hv f 

= E(/I^ElM(^|/)=E(/l P ^vl/)=tr[P,Py] (17) 

/ h v f 

where 

Pv = ElM(V| , (18) 

h v 

is, by definition, the projector onto the Hilbert subspace Hy of confined states 
(i.e. of wavefunctions vanishing out of V). The formula (17) is the starting 
point to carry out a calculation of the MPF at finite volume. The first thing 
to do is to expand (17) as a sum of partition functions at fixed multiplicities, 
i.e.: 

n = E n " (19) 

N 

for a single species gas and: 

n = E Q {Nj} (20) 

{Nj} 

for a multi-species gas, where {Nj} = (N\, . . . , Nk) is a set of particle mul- 
tiplicities for each species j = 1, . . . ,K, defining a channel. Qn or fi{;v-} are 
obtained by summing over all possible values of kinematical variables with 
fixed multiplicities. So, if |/) = \N, {p}) where {p} labels the set of kinemati- 
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cal variables of particles in the state |/), f2jv reads: 

n N = J2( N ,{p}\ p i p v\N,{p}) ■ (2i) 
m 

Likewise, for a multi-species gas, the microcanonical partition function is ex- 
pressed as an expansion over all possible channels: 

= E<TO» Ml^fVITO, M) • (22) 
to 

and ^{Nj} is defined as the microcanonical channel weight. 

The microcanonical channel weights (22) have been calculated in ref. [5] with 
energy-momentum conservation (i.e. using (16) in a multi-particle, first quan- 
tization framework. For a single species ideal spinless gas: 

1 r / n \ n 

QN= m • • • d3pjV 5 \ P ~ £ Pn ) £ II ^(Pp(n) - Pn) (23) 

iV ' J V 71=1 / pe Sjv71=1 

while for a multi-species gas of K spinless bosons: 



{Nj} 



K l Nj 



J=l J Uj=l 



N 



71=1 



j=lpj6 Sat. «j=l 



(24) 



being iV = J2jNj. In eqs. (23), (24) pj labels permutation belonging to the 
permutation group and Fy are Fourier integrals over the region V: 



(25) 



If the volume is large enough so as to allow the approximation: 

F V (P " P') = 7A3 j y d3x e * (P ~ P ' )OC - * S (P ~ P') (26) 



the microcanonical channel weights (24) can be resummed explicitely into the 
microcanonical partition function according to (20) and one obtains [5]: 



lim 



(2tt)< 



j — ic r 

dy° / d 3 y e lP ' y exp 



E 



v 



(27T) 



d 3 p log(l - e- 1 ^)- 1 



(27) 

A full analytical calculation of eq. (27) is possible only for the limiting case 
of vanishing masses (e.g. microcanonical black body). For the massive case, 
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four-dimensional integrations cannot be worked out analytically and one has 
to resort to numerical computation. 

The eq. (27) was implicitely obtained in ref. [4] where the first expression 
of the microcanonical partition function of a multi-species ideal relativistic 
quantum gas was derived as an expansion (20) over channels, by using the 
large-volume approximation (7) from the very beginning. This shows that the 
approximation (26) is indeed equivalent to the (7), as also demonstrated in 
ref. [5]. Noticeably, the MPF definition eq. (14) without any large volume ap- 
proximation involves the appearance of Fourier integrals accounting for Bose- 
Einstein and Fermi-Dirac correlations in the quantum gas, which do not show 
up in the large- volume approximation enforced in ref. [4]. This approach also 
allows to investigate further generalizations when the volume is so small that 
relativistic quantum field effects must be taken into account. 



4 Microcanonical ensemble and field theory 

The calculation of the partition function (14) in a quantum field framework 
brings in new difficulties with respect to the first-quantization scheme. This 
problem has been approached in literature with a functional approach, in- 
spired of the usual grand-canonical thermal field theory [8]. However, these 
calculations aim at the limit of large volumes and are therefore insensitive to 
the difficulties pertaining to the strict requirement of finite volume discussed 
in detail in Sect. 2. As a result, for a free field, the derived expressions are 
equivalent to the formula (27). 

Instead of starting with a functional integration from the very beginning, we 
calculate the microcanonical partition function of a free field by first expanding 
it at fixed multiplicities like in eqs. (19), (21) (or channels, for multi-species 
gas like in eq. (20)), where \N, {p}) are Fock space states with definite particle 
multiplicity and kinematical variables {p}. To carry out this calculation, we 
first need to find an expression of the microcanonical state weight: 

u=(N,{p}\P t P v \N,{p}) . (28) 

By using (16) and choosing \N, {p}) as an eigenstate of the four-momentum 
operator with eigenvalue Pf = J2iPi, the eq. (28) becomes: 

u = 5\P - P f ) (N, {p}\P v \N, {p}) , (29) 

To calculate u> and, by a further integration, Qn we need to know the projector 
Py- Since Py is defined as the projector onto the Hilbert subspace of localized 
states, it can be easily written down in a multi-particle non-relativistic quan- 
tum mechanical (NRQM) framework. As an example, for a non-relativistic 
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spinless single particle, it reads (see also Sect. 2): 



Pv = ElM(M (30) 

k v 

where \ky) is a normalized state of the particle confined in a region V, with a 
corresponding wavefunction ^ v (x) vanishing out of V. The symbol ky stands 
for a set of three numbers labelling the kinematical modes of the confined 
states (e.g. discrete wavevectors, or energy and angular momenta) and the set 
\ky) form a complete set of states for the wavefunctions vanishing out of V. 
The projector (30) can be easily extended to the many-body case and: 

Py= JT \N,{ky})(N,{ky}\ . (31) 
N,{k} 

where the symbol {ky} denotes a multiple set of kinematical modes of the 
confined states while N is the integrated occupation number, i.e. the sum of 
occupation numbers over all single-particle kinematical modes. In the NRQM 
approach these numbers are simply particles multiplicities, implying: 

(N,{p}\N,{k v })^0 iff AT = AT . (32) 

To calculate the microcanonical state weight, thence the MPF, the products 
(N, {ky}\N, {p}) can be worked out on the basis of (32) similarly to what has 
been done in Sect. 2 for a single particle, yielding, for a scalar boson [5]: 

N 

(N,{p}\Py\N,{ P })= ]T II^(Pp(n)-Pn) (33) 

pe Sjv n=l 

where p is a permutation of the integers 1, . . . , N and Fy are Fourier inte- 
grals (25) over the system region V. From the above equation the expression 
of the microcanonical channel weight in eq. (23) follows. We will refer to the 
expression (33) as the NRQM one, meaning that is has been obtained in this 
first-quantized multi-particle NRQM framework, where, for instance, particles 
and antiparticles are simply considered as different species and their contri- 
butions factorize. 

One could envisage that a projector like (31), written in terms of Fock space 
states, could be simply carried over to the relativistic quantum field case, 
where \N, {ky}) are states of the localized problem, obtained by solving the 
free field equations in a box with suitable boundary conditions. Yet, some 
difficulties soon arise. First of all, whilst in NRQM the single-particle local- 
ized wavefunction \ky) and the free plane wave state |p) live in the same 
Hilbert space, in quantum field theory the localized and the non-localized 
problem are associated with distinct Hilbert spaces. Thus, unlike in NRQM, 
it is not clear how to calculate a product like (N, {k v }\N, {p}). Secondly, even 
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if there was a definite prescription for it, it should be expected that the inte- 
grated occupation numbers of the localized problem do not coincide with actual 
particle multiplicities unless the volume is infinite. To understand this point, 
one should keep in mind that properly called particles arise from solutions of 
the free field equations over the whole space and that the hamiltonian eigen- 
states of the localized problem are conceptually different. Consequently, the 
integrated particle number operators in the whole space should differ from in- 
tegrated number operators within the finite region. Hence, unlike in NRQM, 
a state with definite integrated occupation numbers N (we purposely refrain 
from calling them particle numbers) should be expected to have non vanish- 
ing components on all free states with different numbers of actual particles, 
namely: 

\N)v = a o x\0) + a lt fi\l) + ... + a^\N) + ... (34) 

where a. ^ are non-trivial complex coefficients, and (32) no longer holds. Only 
in the large volume limit one expects that the integrated occupation numbers 
coincide with actual multiplicities and eq. (32) applies. This kind of effect is 
pointed out in the introduction of Landau's book on quantum field theory 
[7]: when trying to localize an electron, an electron-positron pair unavoidably 
appears, meaning that the localized single "particle" is indeed a superposition 
of many true, asymptotic particle states. Another relevant manifestation of 
this difference which is probably more familiar, is the Casimir effect, which is 
related to the difference between the true vacuum state |0) and the localized 
vacuum state \0)v- Hence, all formulae derived under the approximation (32) 
are asymptotic ones, valid in the limit V — > oo but not at strictly finite volume. 
Thus, one should expect significant finite-volume corrections to the eqs. (33) 
and the ensuing MPF (23) in a quantum field treatment. 

If we want to give an expression like (N, {ky}\N, {p}) a precise meaning in a 
quantum field framework, we first need to map the Hilbert space Hy of the 
localized problem into the Hilbert space H of the free field over the whole 
space. This can be done by mapping the field eigenstates and operators of H 
into Hy in a natural way as: 



#(x)^ v — > (x) H 

|V(x))h v — (35) 

This allows writing linear, non-bijective, Bogoliubov relations expressing the 
annihilation and creation operators of the finite region problem as a function 
of those of the whole (real scalar) field (see Appendix A for derivation): 

* = I d ' p F{k - P) ^0 P tt * + F(k - -">^ < (36) 

where k are triplets of numbers labelling kinematical modes, just like the 
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aforementioned ky, £k is the associated energy, e p = VP 2 + m ' 2 an d: 

f(k,p) = (2^{ d3x<(x)e,P ' X (37) 

Uk being a complete set of orthonormal wavefunctions for the finite region. 
A remarkable feature of relativistic quantum fields is that, unlike in NRQM, 
the localized annihilation operators have non-vanishing components onto the 
creation operators in the whole space, as shown by (36). This confirms our ex- 
pectation that a localized state with definite integrated occupation numbers is 
a non-trivial linear combination of states with different particle multiplicities. 
Expectedly, as the volume increases, these components become smaller and in 
the infinite volume limit one recovers = a p (see Appendix A). 

Starting from the Bogoliubov relations (36), it should be possible, in principle, 
to calculate the coefficients in eq. (34), thence the microcanonical state weight 
(29) by using the expansion (31). In fact, we do not really need to do that. 
It is more advantageous, as pointed out in ref. [6], to write the projector Py 
in terms of field states rather than occupation numbers of field modes within 
the finite region. Indeed, in the general definition in eq. (17): 

Pv = ElM(M (38) 

hy 

the states \h v ) are a complete set of states of the Hilbert space of the localized 
problem Hy, where the degrees of freedom are values of the field in each 
point of the region V, i.e. {^(x)} | x e V. Therefore, the above projector is a 
resolution of the identity of the localized problem and can be written as (for 
a real scalar field): 

P v = f Zty#>M (39) 

where = Cg> x |-(/>(x)) and T>ip is the functional measure; the index V means 
that the functional integration must be performed over the field degrees of 
freedom in the region V, that is T)i\) = Ilxev c ^/ ; ( x )- The normalization of 
the states is chosen so as to (^(x)|^//(x)) = <5(t/>(x) — tp'^x)) to ensure the 
idempotency of Py. If we now want to give expressions like: 

({NjyApMPvmhip}) (40) 

a clear meaning, we should find a way of completing the tensor product in the 
projector (39) with the field states outside the region V such a way the scalar 
product can be performed unambiguously. 

Unfortunately the answer to this question is not unique and the projector 
can be extended to H in infinitely many ways. What is important is that the 
result of the calculation is independent of how the projector has been extended. 
Thus, at the end of the calculation, one has to check whether spurious terms 
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appear, possibly divergent, depending explicitely on the chosen state of the 
field outside V and these terms must be subtracted away. In general, all terms 
depending on the degrees of freedom of the field out of V must be dropped 
from the final result. 

In this work, we will extend the projector with eigenstates of the field, where 
the field function i/j(x) is some arbitrary function outside the region V. Thus, 
the projector Py (39) is mapped to: 

P v = f VtP\i/j){iP\ |V) = ®* 6 v|V(x)) |V(x)) (41) 

J v 

where the index V still implies T>%p = Ilxev dip(x). We will see that, with the 
definition (41), spurious terms depending on the degrees of freedom outside 
V do arise indeed, but that they can be subtracted "by hand" in a consistent 
way. 



5 Single particle channel 

We will start calculating the expectation value (p\Pv\p) of a single parti- 
cle channel in the simple cases of neutral and charged scalar fields. This is 
preparatory to the general case of multiparticle states in Sect. 6. 

5. 1 Neutral scalar field 

We consider a gas made of one type of spinless boson described by the free 
real scalar fielqjj (in Schrodinger representation): 

¥(x) = / ^= ( a(p) e*« +a\p) e^ x ) (42) 

(27r)2 J y/2e v ' 

where e = Vp 2 + rri 2 is the energy, p is the modulus of the three-momentum 
and the normalization has been chosen so as to have the following commutation 
rule between annihilation and creation operators: 

[a(p),a\p')} = S\p-p'). (43) 

We start writing the one-particle Fock state \p) in terms of creation and an- 
nihilation operators acting on the vacuum: 

( P |Pvb> = <0|«(p)Pvo t (p)|0>. (44) 

3 Henceforth, the capital letter ^ will denote field operators while for field functions 
we will use the small letter if;. 
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Since Py is denned, according to (41) as a functional integral of eigenvectors 
of the field operator it is convenient to express creation and annihilation 
operators in terms of the field operators. We shall use following expressions 
which are the most appropriate for our task: 



at 



<0|a(p) = <0| ~~~~tt / d 3 x e-^ x v / 2i^(x) 
(27r) 2 J 

(P)|0> = 7^7T / d 3 x e«P*v^*(x) |0) . (45) 

(27TJ 2 J 



These expressions can be easily checked by plugging, on the right hand side, 
the field operators in (42). By using eq. (45) in (44): 

(0\a(p)P v a^(p)\0} = J d 3 x J d 3 x' e tp ^-^2e (0|*(x')Py#(x)|0) (46) 

It can be easily checked now that the rightmost factor in the above equation 
turns out to be (by using the definition (41)): 

(0|*(x')Py*(x)|0) = / 1(^10)1 WM*) • (47) 

where ip(x) and ?/>(x') are field functions or the eigenvalues of the field operator 
relevant to the state that is: 

*(x)|^) = |^(x) . (48) 

It is possible to find a solution of the functional integral (47) by first consid- 
ering the infinite volume limit, when the projector Py reduces to the identity. 
In this limiting case, the functional integral in (41) is now performed over 
all possible field functions and eq. (41) becomes a resolution of the identity; 
(0|\I / (x / ) , I'(x)|0) is just the two-point correlation function that we write, ac- 
cording to (47): 

(0|*(x')*(x)|0> = Jvij; |(^|0)| 2 ^(x>(x) (49) 

The product (-0 j 0) is known as the vacuum functional and reads [12], for a 
scalar neutral field: 

^|0)=Arexpj-i|d 3 x 1 |d 3 x 2 ^(x 1 )ii(x 1 -x 2 )^(x 2 )} (50) 

where M is a field-independent normalization factor, which is irrelevant for 
our purposes. The function K(x' — x) is called kernel and fulfills the equation 
[12]: 

J d 3 x' e- ip ' x ' K(jd - x) = 2e e~ ip x (51) 
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whose solution is: 



K(x' - x) = J d 3 p e~ ip ' (x '- x) 2e . (52) 

The functional integral (49) is therefore a gaussian integral and can be solved 
by using the known formulae for multiple gaussian integrals of real vari- 
ables [12]: 



2N r l r r 

I 2N = JV^J II^W ~2J d " Xl j /d 3 x 2 ^(x 1 )K(x 1 -x 2 )^(x 2 ) (53) 

i=l 

= ^2 H K~ l (paired variables) 



pairings pairs 
Of t,l,---,t,2N 



where paired variables means couples whose difference — £j (or, what 

is the same, £j — £j as K" 1 is symmetric) is the argument of K~ l . The factor 
J is just the normalization of the vacuum state I = (0 1 0) which is set to 1. 
The inverse kernel K~ x can be found from its definition: 



d 3 x / J ^(y-x / )^ 1 (x / -x) = 5 3 (x-y) Vx,y (54) 
leading to: 

K~\* - x) = JL J ^ e-^*'-*) = (0|*(x')*(x)|0> (55) 

The last equality comes from (53) and (49) in the special case N = 2 or can 
be proved directly from field Fourier expansion (42). 

We are now in a position to solve the functional integral (47) at finite volume. 
First, the functional integration variables are separated from those which are 
not integrated, i.e. the field values out of V (V denotes the complementary of 
V): 



| M0> | 2 = \N\ 2 exp -- J d 3 Xl J d 3 x 2 V(x!)X( Xl - x 2 )^(x 2 ) 
= |A/T exp [-1 J d 3 x x J d 3 x 2 V(xi)X( Xl - x 2 )^(x 2 ) 



x exp 
x exp 



-- J d 3 Xl J d 3 x 2 V(x 1 )X(x 1 - x 2 )^(x 2 ) 
- |^d 3 Xl |_d 3 x 2 ^(x 1 )K(x 1 -x 2 )V(x 2 



(56) 



where we have taken advantage of the symmetry of the kernel K. The (56) 
is a gaussian functional, with a general quadratic form in the field values in 
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the region V; it can be integrated in (47) according to standard rules [12], 
yielding: 



r^i- 1/2 

exp 



2vr 



/ |<0|^>|WMx) = ^ 1 (x / ,x)|AT| 2 det 
X - J d 3 Xl d 3 xi J d 3 x 2 d 3 x' 2 K v \^ x / 1 )^(x 1 - x 2 )K(x; - x' 2 )V(x 2 )V(x 2 ) 
-^/_d 3 Xl |_d 3 x 2 ^(x 1 )^(x 1 -x 2 )^(x 2 ) 

= ^(x',x) / |(0|^)| 2 = K 1 7 1 (x',x)(0|P y |0) (57) 

The function Ky 1 is the inverse of K over the region V, namely the inverse 
of: 

K(x! - x)0 y (x / )e y (x) = K v (v! - x) , (58) 
the function 0y(x) being the Heaviside function: 



,1 if x G V 

Ov(x) = <( (59) 
otherwise. 

The inverse kernel .fTy 1 fulfills, by definition, the integral equation: 

f d 3 x' K v (y - x')K v \x', x) = 5 3 (x - y) VxjGV (60) 
«/ v r 

Note that, because of the finite region of integration, the inverse kernel may 
now depend on both the space variables instead of just their difference. Also 
note that Ky 1 is real and symmetric, being Ky real and symmetric. Therefore, 
the result of the functional integration yields the simple formula: 

(0|*(x')Py*(x)|0> = ^ 1 (x',x)(0|P y |0) (61) 

where the factor (0|Py|0) is a positive constant which we will leave unex- 
panded. 

Altogether, the presence of the projector Py in the eq. (61) modifies the two- 
point correlation function by introducing a constant factor (0|Py|0) and re- 
placing the inverse kernel K^ 1 with a different one Ky 1 . It can be easily 
proved, by using the general formulae of gaussian integrals, that this holds 
true in the more general case of many-points correlation function. In fact, the 
(53) holds for general quadratic forms in the field ip and so the eq. (61) can 
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be generalized to: 

N 2N 

(0|n*K)Py II *W|0> = <0|Pv]0> ]T I] ^'(Paired var.) . 

n=l n=jV+l pairings pairs 

of Xl,...,X2N 

(62) 

The remaining task is to calculate the inverse kernel Ky 1 by means of (60). 
In fact, we will look for a solution of the more general equation: 

f dV K(y - x')K v \x', x) = 5 3 (x - y) V x e V, y (63) 

with unbounded y. It is clear that a solution Ky 1 of equation (63) is also 
a solution of (60) because K v equals K when y £ V. Setting y unbounded 
allows us to find an implicit form for Ky 1 . In fact (63) implies: 

1 r , e ipx 

w( dV{ (64) 

which is obtained multiplying both sides of (63) by e* py / (27r) 3 and integrating 
over the whole space in d 3 y. 

We are now in a position to accomplish our task of calculating (p|Py|p). By 
plugging (61) into (46) we get: 

(0|o(p)Pvat(p)|0> = ^ J d 3 x J dV e^( x ~ x ') 2e Ky 1 ^', x)(0|P v |0> . 

(65) 

The integration domain in (65) can be split into the region V and the com- 
plementary V for both variables. The inverse kernel Ky 1 is not defined out of 

V and can then be set to an arbitrary value, e.g. zero. Otherwise, even if one 
chose a non vanishing prolongation of Ky 1 , an integration outside the domain 

V would involve the degrees of freedom of the field out of V and, according to 
the general discussion at the end of Sect. 4, the contributing term should be 
dropped. Therefore, retaining only the physically meaningful term, the (65) 
turns into: 

(0|a(p)P y at(p)|0) = j v d 3 x |^dV e^ x - x ') 2e Ky\x!, x)(0|P v |0) . 

(66) 

In the above equation one can easily recognize the complex conjugate of the 
left hand side of (64). Hence, replacing it with the complex conjugate of the 
right hand side, one gets: 

(0|a(p)Pvat(p)|0> = j^- 3 /^d 3 x(0|P v |0) = ^^{0|P V |0> (67) 
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which is the same result of NRQM in (33) in the simple case N = 1, times 
a factor (0|Py|0). This factor still contains a dependence on the field degrees 
of freedom out of V, according to the projector expression in eq. (41), which 
should disappear at some point. However, we will see that this factor appears 
at any multiplicity and therefore becomes irrelevant for the calculations of the 
statistical averages. 



5.2 Charged scalar field 



The calculation done for a neutral scalar field can be easily extended to a 
charged scalar field. The 2-component charged scalar field in Schrodinger rep- 
resentation reads: 



¥(x) = -U / ^ ( a(p) +&t(p) e -^ x ) (« 
(27r)2 J V 2s v ' 

*t( x ) = -L r [<PP( b(p) e^ x +at(p) e"^ 



where a, and 6, 6^ are annihilation and creation operators of particles and 
antiparticles respectively. They satisfy commutation relations: 



Hp), a\p')\ = [b(p), b\p')\ = 5 3 (p - p') (69) 
[a(p),b(p')] = [a\p),b(p')]=0 

Likewise, the fields obey the commutation relations: 

[*(x),*t(y)] = . (70) 

and it is then possible to construct field states \tp,ijr). The projector Py can 
be written as: 

P V EE / P(^)|^t)(^t| (71) 

with suitable state normalization and arbitrary field functions if>(x) out of the 
region Similarly to eq. (45), one can write: 



4 The functional measure in the equation (71) reads Y\ x dip(x)dip* (pcj/in. Anyhow, 
its explicit form is not important for our purposes. 
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(0\a(p) = (0| [ d 3 x e" ip - x Vte *(x) (72) 

(27TJ2 •/ 

0+07)10) = — ^— 3 / d 3 x e^v^e^^x) |0) 

(27f)2 J 

(0|6(p) = <0| -i-r / d 3 x e^v^^x) 
6t( p )| > = — L_ /" d 3 x e i P' x v / 2i^(x) |0) 

(27r)2 J 

The chain of arguments of the previous subsection can be repeated and the 
functional integral: 

N N N 

n=l n=l JV n=l 

(73) 

found to be a multiple gaussian integral. Letting p be a permutation of the 
integers 1, . . . , N, the integration on the right hand side of eq. (73) yields: 

N N N 

<0| II V^PV II = (0|Pv|0> E II ^(^(n)) (74) 

n=l ra=l pGSjv n=l 

which differs from the corresponding expression for the real scalar field because 
now the field is complex and ip can only be coupled to ip^ [13]: 



N r 

[ V(^,4>) J] (V(^n)V' t (0)exp - / d 3 Xl / d 3 x 2 ^( Xl )^( Xl _ X2)V; ( X2 ) 

J n=l 1 J J 

N 

= / o E II " (75) 

However, the functional integral involving only two fields ip and ^ yields the 
same result as for neutral particles. Thus, the kernel K is still the same and 
so is the integral equation (64) defining Ky 1 . The expectation value of Py 
in a state with only one particle (or antiparticle) will also be the same as in 
eq. (67), that is: 

(0\a(p)P v a\p)\0) = (0\b(p)P v b\p)\0) = J^(0|P y |0> . (76) 



6 Multiparticle channels 



We have seen in the previous section that the expectation value (p\Pv\p) f° r a 
single spinless particle is the same obtained in a NRQM approach [5] times an 
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overall immaterial factor (0|Py|0). In this section, we will tackle the calcula- 
tion of the general multiparticle state. We will see that, by using the projector 
definitions in eqs. (41), (71) and subtracting the spurious contributions stem- 
ming from external field degrees of freedom, the final result is still the same 
as in the NRQM calculation times the factor (0|Py|0). 



We will first address the case of N charged particles. 



6. 1 Identical charged particles 



We will consider a state with N identical charged particles; for N antiparticles 
the result is trivially the same. 

We want to calculate: 



N 



N 



(N, {p}\P v \N, {p}) = (0|n«(Pn)PyIl at (Pn)|0) . 



(77) 



Since: 



[a\ tft] = [ a> = [ 6> ^t] = [ 6 t > ^r] = o 



(78) 



one can replace creation and annihilation operators with their expressions in 
(72) disregarding the position of the operators with respect to the vacuum 
state. In formula: 



N 



(n, { P }\p v \n, M)=n 



71=1 



^- J d 3 x n J d 3 x'„ e-^-^-O 2e Ti 



(2nY 

j 



n=l 



n=l 



N 

n 

n=l 



(27T) 
N 



J d 3 x n J d 3 x'„ e-^-^-O 2e r , 



x E II^v^^XOlPvlO), 

pSSjv n=l 



(79) 



(80) 



where is the permutation group of rank N. Now, like for the single particle 
case, we restrict the integration domain to V in (79) in order to get rid of 
external degrees of freedom and, by repeatedly using eq. (64), we are left 
with: 
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(N,{p}\Py\N,{ P }) _ " 

(OlPvlO) ~ tkjA 

N 

= e n 

pSSjv n=l 



Hence, using the definition (25): 



1 / d 3 x' e~ ip "- (x p(") _<) 
(2tt)3 Jv n 



e -*(Pp-l(„)-Pn)-x; 



iV 



(iV, { P }\P V \N, { P }) = E II ^(Pp(n) - P„)(0|P v |0) 

peS N n=l 



(81) 



(82) 



which is exactly the expression (33) obtained in NRQM for N identical bosons, 
times the factor (0|Py|0). 



6.2 Identical neutral particles 



The case of N identical neutral particles is more complicated because of the 
possibility of particle pairs creation. This is reflected in the formalism in the 
occurrence of many additional terms in working out the eq. (77). We start 
with the calculation for a state with two neutral particles of four-momenta p\ 
and pi respectively, generalizing to N particles thereafter. In terms of creation 
and annihilation operators: 

(pi,p 2 |Pybi,P2> = (0\a( Pl )a(p 2 ) P v a* (p 2 )a* ( Pl )\0) (83) 

which we can rewrite using eq. (45) for a{pi) and a\pi) as: 



<Pi,P2|Pvbi,P2> = ^3 / d 3 x x J d'M e-**"**™') 

x 2e 1 (0|*(x 1 )a(p 2 ) P v ^(^(x'JlO) . (84) 

In order to use the expression (45) for annihilation and creation a(p2) and 
a^(p2) operators we have to get them acting on the vacuum state, hence they 
should be moved from their position in eq. (84) outwards. This can be done 
by taking advantage of the following commutation rules: 



[*(x),a(p)] 
[a(p)t,¥(x); 



Q—ipx 



(27r)2V2e 

gip-x 



(27T)5v^i 



(85) 
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Using (85) and then plugging (45) in eq. (84) we get: 



(Pi,P2|Py|Pi,P2) = 

d 3 Xl f d 3 x; e- ipi ( Xl - Xl ')e- tp2 -( x ^ x 'i)(0|Py|0) 



e 2 (27r)6 
2ei 



J d 3 Xi y d 3 x; /" d 3 x 2 e - ipi -( Xl - X1 ') (0|*(x 2 )P y *(x 1 )|0) 
f d 3 Xi y dV^ d 3 x' 2 e- ipi '( Xl - X1 ') e- ip2 < Xl ~ x '^ (0|^(x , 1 )P y ^(x / 2 )|0) 



(2k) 
2ei 



4£i£2 

+ (2tt) 6 

x(0|^(x 2 )^(x 1 )P y ^(x / 1 )^(x / 2 )|0) 



y d 3 Xl y d 3 x; y d 3 X2 y dv. 



-ipi-(xi-xi') „-ip2-(x2-x 2 ') 



(86) 



where advantage has been taken of the fact that [Py, ^(x)] = 0. The eq. (86) 
has four different terms, among which only the last was present in the case 
of 2 charged particles in Subsect. 6.1. Instead, the first three terms arise from 
contractions of the annihilation and creation operators with field operators on 
the same side with respect to Py. By using (62) and replacing all unbounded 
integrations with integrations over V to eliminate spurious degrees of freedom 
of the field, the eq. (86) yields: 



(Pl,P 2 |Py|Pl,P 2 ) _ £l f , 3 f ,3 / p-ipi-Cxt-X!') e -<P2-(xi-xl) 

(0|Py|0) ~ e^nyJv^^Jv 1 

j^d 3 xi j^d^a j^dVj e - ipi '( Xl - X1 ') e - ip2 '( X2 - x 'i) K l ; 1 (x 2 ,x 1 ) 

(2tt) 6 

, 4£i£ 2 



J d 3 Xl J dM J d 3 X 2 e -iPii*i-xi') e -*P2-(xi-x'2) Ky\^*! 2 

I d 3 Xl / dV, / d 3 x 2 / d 3 x' e - ipi -( Xl - X1 ') e - tp2 ( X2 - X2 ') 
Jv Jv Jv Jv 



X 



Ky^^Ky 1 ^, X 2 ) + Ky 1 ^, X^ICy^, x'J 



+iv^ 1 (x 1 ,x 2 )K v ; 1 (x 2 ,x / 1 ) 



(87) 



We can now use (64) to integrate away the inverse kernel in eq. (87). For 
the second term on the right hand side, we choose to integrate exp(— ip 2 • 
x 2 )i^y 1 (x 2 , xi) in x 2 and for the third term exp(ipi • x / 1 )fCy 1 (x / 1 , x' 2 ) in x[. 
At the same time, we perform the integration of the last term in an arbitrary 
couple of variables with indices 1 and 2, obtaining: 
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(0|P y |0) (2ir) G Jv Jv 1 e 2 

1 /■ /• e -ip2-(xi-x' 1 ) 

/ d 3 x x / d 3 x' e-'P 1 ^ 1 -^') 2ei 

(2n) 6 Jv Jv 1 2e 2 

_ f H 3 ^ / H 3 y' p-«Pl'( x l- x 2') p -iP2-(xi-x 2 ) 

{2^Jv 1 Jv 2 



V 2 , 

-j(pi-p 2 )-(xi-x 2 ) 



+— f 77~~~T7r / d 3 ^! / d 3 x 2 e 

(27r) 6 (2tt) 6 Jv Jv 



+ 



^— / d 3 Xl / d 3 x' e - ipi '( Xl - X2 ') e - ip2 '( Xl - x 2) (88) 
l7r) b Jv Jv 



(2tt)< 

Four terms in the above sum cancel out and we are left with: 

(pi,P2\Py\pi,P2) V 2 , I p / ^ v ,2 V- A F /„ \ / 8Q x 

<0|P V |0> = Wf + 1 FWPl - P2) 1 = S 2 R MPp(n) " { ) 

This is the same result one would obtain for two identical bosons in a NRQM 
framework [5]. The second term on the right hand side accounts for the well 
known phenomenon of Bose-Einstein correlation in the emission of identical 
boson pairs. 

Looking back at the whole derivation, we find that the terms arising from 
pairings of field variables on the same side with respect to Py have cancelled 
with the terms stemming from commutation of annihilation and creation op- 
erators with field operators; the only surviving terms are N\ pairings of field 
variables on different sides of Py, just like in the case of charged particles. This 
cancellation property holds and so, the formula (82) applies to the case of ./V 
neutral particles as well. A proof based on the form of the thermodynamic 
limit V — > oo is given in Appendix B. 



6.3 Particle- antiparticle case 



For a state with one particle and one antiparticle, the expectation value of Py 
reads: 

(0\a( Pl )b(p 2 ) P v tf( P2 W( Pl )\0) (90) 
implying, in view of (72): 

(0\a( Pl )b(p 2 ) Py tf(p 2 )a\ Pl )\0) 

= J d 3 x x J d 3 x' x e-^ ( Xl - Xl ')2£ 1 (0|M/(x 1 )6b 2 )Py 6 t (p 2 )^ t (x' 1 )|0)(fll) 

Like for neutral particles, the b and operators are moved outwards to get 
them acting on the vacuum state by using the commutators: 
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[tt(x),6(p)] 



[6(p)t,^( x )] = e 7 ( 92 ) 
(27r) 2 V2e 

By using the eq. (72) and, like in the previous subsection, restricting the 
integration to the region V to eliminate external degrees of freedom, we obtain: 



(0\a( Pl )b(p 2 ) P v b\p 2 )a\ Pl )\0) 

■ J d 3 *! J d 3 ^ e- ipi -( Xl - Xl ')e- ip2 -( Xl - x 'i)(0|P y |0) 



e 2 (2ir 



6 

2ei 



f d 3 Xl / d 3 x; / d 3 x 2 e - ipi '( xl - X1 ') e- ip2 '( X2 - x i) (0|^ t (x 2 )P y ^(x 1 )|0) 
(2tt) 6 X ^ X ^ X ^ e_<Pr(xi_Xl ' ) e- ip2 '( Xl - x 2) (0|^ t (x' 1 )P y ^(x , 2 )|0) 



(2k) 
2ei 



4eie 2 



• •• : / d 3 Xl / d 3 x' x / d 3 x 2 / d 3 x 2 e ^ 1 <> 

Jy Jv 



jpi-(xi-xi') „-ip2-(x2-x 2 ') 



(2^7^ Vy 2 Jv 2 

x (0|*t(x 2 )tt( Xl ) P y M/t(x / 1 )M/(x' 2 )|0) (93) 

whose result is (87) except the (last-l)th term, because now: 



(0|M/t(x 2 )M/( Xl ) P v M/t( x ;)M/(x 2/l , 

(0|P V |0> (94j 

which is a consequence of the general expression (74). Again, four out of five 
terms cancel out in eq. (93) and we get: 

(0\a( Pl )b(p 2 )P v tf(p 2 )al( Pl )\0) V 2 



(0\P V \0) (2vr)6 

where the missing term with respect to the neutral particle case (89) is the 
one involving permutations; this is a natural result because particles and an- 
tiparticles are, of course, not identical. Hence, the expectation value of Py on 
a particle-antiparticle pair shows a remarkable factorization property, that is: 

(0\a( Pl )b(p 2 ) P v tf(p 2 W( Pl )\0) = (0\a(p 1 )Pva^p 1 )\0){0\b(p 2 )P v b^p 2 )\0) 

(96) 

where the = sign holds provided that the integrations are restricted to the 
region V. 

Extrapolating to the most general case of N + particles and iV_ antiparticles, 
it can be argued by using the limit V — > oo that the factorization of the 



24 



microcanoncal state weight holds (see Appendix B) at any multiplicity, i.e. 
particles and antiparticles behave like two different species: 



({N + }, { P+ }, {N_}, {p_}\P v \{N + }, { P+ }, {AL}, {p_}) 

N + N— 

: E II MP P+ (n + ) - Pn + ) E II ^(Pp-(n_)-Pn_)(0|Py{OX) 

p+eSiv, n+=l p-eSjv_ n_=l 



7 Summary and discussion 



On the basis of eqs. (19), (21), (16) and (82), which applies to charged as 
well as neutral particles, and taking into account that the states \N, {p}) have 
been chosen to be eigenvectors of four-momentum, we can finally write down 
the full expression of the microcanonical partition function of a relativistic 
quantum gas of neutral spinless bosons as: 



OO 1 

^=(0|Py|0)E 



AT=0 



N 

n / d3 pn 

71=1 J 



N 



n=l 



N 

x E II / d 3 xexp[ix- (p p(n) -p n )] 

pGSjv n=l JV 



(98) 



with P = (M, 0) and the factor 1/N\ has been introduced in order to avoid 
multiple counting when integrating over particle momenta. Similarly, on the 
basis of eqs. (20) (22) and (97), the microcanonical partition function of a 
relativistic quantum gas of charged spinless bosons can be written as: 



n = <o|Pv|o> £ 



N + ,N-=0 + 



N++N- 



n " / d 3 P, 



N++N- 
S 4 \P- E Pn 



n=l 



x E II / d 3 xexp[ix- (p p+(n+ ) - p n+ )] 

p+eS N+ n+=l v 
N ~ r 

Ell / d 3 X6Xp[ix- (p p _(n_) -Pn_)] • 



P_GSat_ 



(99) 



The generalization to a multi-species gas of spinless bosons is then easily 
achieved: 
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n = <o|p v |o> £ 

{Nj} 

k N. 



K 1 



N 



n=l 



s 4 p-Ep, 

\ n- 

x fl E II L d 3 xexp[^x • (p p . (n .) - p n ,)] . 

j=l pi& S N , n,=l 



(100) 



where iV = The formulae (98), (99) and (100) are our final result. 

The finite volume Fourier integrals in the above expressions nicely account 
for quantum statistics correlations known as Bose-Einstein and Fermi-Dirac 
correlations. We stress once more that for a charged quantum gas, particles 
and antiparticles can be handled as belonging to distinct species and they 
correspond to different labels j in the multi-species generalization of (100). 



The expression of the MPF (98) is the same as obtained in a NRQM calculation 
in ref. [5], quoted in this work in (24), times an overall factor (0|Py|0) which is 
immaterial for the calculation of statistical averages in the microcanonical en- 
semble. More specifically, the expectation value of Py on a free multi-particle 
state (see eqs. (82) and (97)) is the same as in the NRQM calculation (33) 
times (0|Py|0). This result has been achieved enforcing a subtraction prescrip- 
tion, namely that all terms depending on the degrees of freedom outside the 
system region V must be subtracted "by hand" in all terms at fixed multi- 
plicities. The factor (0|Py|0) is still dependent on those spurious degrees of 
freedom, according to the Py definition in eq. (41), but this does not affect 
any statistical average because it always cancels out. In the thermodynamic 
limit, this factor tends to 1 as Py — > I and the large volume limit result known 
in literature [4] (i.e. eq. (27)) is recovered. 



This result looks surprising in a sense because one would have expected, a pri- 
ori that quantum relativistic effects would affect the statistical averages bring- 
ing in a dependence on the ratio between the Compton wavelength and the 
linear size of the region, getting negligible at small values, i.e. when V — > oo. 
This is because the condition (32) from which the MPF expressions (23,24) 
ensue, no longer holds in quantum field theory at finite volume and only ap- 
plies with good approximation at very large volumes, as discussed in Sect. 4. 
However, in the calculation of statistical averages, only the projector Py enters 
and this implies^the summation over all states at different integrated occupa- 
tion numbers \N)v, according to (31). So, even though the coefficients in the 
expansion (34) are different from zero and do depend on the volume, it turns 
out that summing all of them at fixed N, one gets the same result as in the 
NRQM approximation (32). In formula, by using (31): 
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(N,{p}\Pv\N{p})= £ \(N,{p}\N,{k v }) 



N,{k v } 



E I(o|^,{M>I 2 EK^MI^{M>Ia 



(101) 



N,{k v } i k v} 



where the index NR means that one should make a non-relativistic quantum 
mechanical calculation. 

This calculation can be extended to the most general case of the microcanon- 
ical ensemble of an ideal relativistic quantum gas by fixing the maximal set 
of observables of the Poincare algebra, i.e. spin, third component and par- 
ity besides the energy-momentum four-vector. This will be the subject of a 
forthcoming publication. 
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A Bogoliubov relations for a real scalar field 

They can be derived by first expressing the localized annihilation operator 
as a function of the field ^ and its conjugated moment For the localized 
problem we have: 



where k is a vector of three numbers labelling the modes in the region V, e k 
is the associated energy and is a complete set of orthogonal wavefunctions 
over the region V: 





XX( x W(x') = 5 3 (x-x) 



k 




(A.2) 



and vanishing out of V . Inverting the (A.l) we have: 




(A.3) 
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which are valid at any time. We now enforce the mapping (35) and replace 
the localized field operators at t — 0, i.e. in the Schrodinger representation, 
with those in the full Hilbert space. In other words, we replace in (A. 3): 



* (x) ^(2^/ d3pZ vr e " p ' Xa > +c - c - (A - 4) 



where e p = Vp 2 + m 2 , obtaining: 



where: 

F(k, p) = j^j- 2 jf d 3 x<(x) e^ x (A.6) 

that is the Bogoliubov relations (36). We observe that the localized annihila- 
tion operator is a non-trivial combination of annihilation and creation opera- 
tors in the whole space. However, the term in (A. 5) depending on the creation 
operator a* vanishes in the thermodynamic limit, as expected. This is most 
easily shown for a rectangular region, where one has «k( x ) — exp[ik • x\/y/V ; 
k is given by the eq. (9) and £ k = y/k 2 + m 2 . Hence: 

lim F(k, -p) oc lim / d 3 x e -^ p+k > x oc 5 3 {p + k) (A.7) 

and, consequently e^— e p — > 0, which make the second term in (A. 5) vanishing. 

We can use the eq. (A. 5) to work out linear relations between Fock space 
states of the localized problem and asymptotic Fock space states. These can 
be obtained enforcing the destruction of the localized vacuum state: 

a k \0) v = 

and writing: 

|0)v = ao|0) + /d 3 pai(p)|p) + j d 3 Pl d 3 p 2 a 2 (pi, p 2 )|pi, p 2 > + • • • (A.8) 
Working out such relations is beyond the scope of this paper. 
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B Multiparticle state with N particles 



In order to prove eq. (82) for neutral particles, the first step is to realize that 
a general expansion of: 

(0\a l ...a N P v a t N ...a\\0) (B.l) 

(dj is a shorthand for a{pif) must yield, on the basis of what shown for the 
case of 2 particles, a sum of terms like these: 

/ d 3 Xl .../ d 3 x N n e ±i ^f(e 1 ,...,e N )(0\P v \0) (B.2) 

^ ^ n,m=l 

where / is a generic function involving a sum of ratios or product of any 
number of particle energies and ± stands for either a + or a — sign. That (B.2) 
ought to be the final expression can be envisaged on the basis of a repeated 
application of the equations (45), (85), (62), (64) in turn. The eq. (B.2) can also 
be written as: 

J d 3 Xl ...J d 3 x N e^>"- Xl . . . e 1 f{e u e N )(0\P v \0) (B.3) 

where J2' stands for an algebraic sum with terms having either sign. 

If we now take the thermodynamic limit V — > oo of (B.l), one has Py — > I, 
thence: 



lim (0|ai . . . aArPyajy . . .a\\0) = (0|ai . . . a^a^ . . . a\\0) 

= E Il^(Pn-Pp(n)) (B.4) 

pGSjv n=l 

This tells us that the function /'s in each term (B.2) must reduce to a triv- 
ial factor 1 because they would otherwise survive in the thermodynamic limit, 
being a factor depending only on the p n 's. Moreover, since the thermodynamic 
limit involves only Dirac's 5's with differences of two momenta as argument, 
there can be only difference of two momenta as argument of exponential func- 
tions in (B.3). Finally, by comparing the eq. (B.3) with the (B.4), we conclude 
that the only possible expression at finite V is: 

(0|a 1 ...a^aj v ...a t 1 |0) = (0|P^|0> ]T f[ — L- / d 3 x n e^-^))'*" 

pe s N n=i y Z7r ) Jv 

(B.5) 

which is precisely (82). 
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For a state with iV + particles and AL antiparticles, the validity of eq. (97) can 
be argued with a similar argument, i.e. by constraining the form of general 
terms like (B.2) taking advantage of the limit V — > oo. In the case of particles 
and antiparticles, the thermodynamic limit tells us that: 

lim (0|ai . . . a N bi . . . b N _P v a i N . . . a\tf N _ . . . b\\0) 

= (0|ai . . . a N+ h . . . b N _a ] N+ . . . a[V N _ . . . 6[|0) 

n + v_ 

= E II 53 (Pn+ - Pp+(n+)) E II ^(Pn- ~ Pp-(«-)) ( B ' 6 ) 

p+e5jv + n+=l p-eSjv_ n_=i 

and this determines the form of the expression: 

(0|ai . . . a N+ bi . . . b N _P v a) N+ . . . a[V N _ . . . &j|0) 

at finite V to be (97). The absence of integrals mixing particles with antipar- 
ticles momenta such as: 

J d 3 Xi e l( - Pn +~ Pn -> Xi 
is owing to the absence of such differences as arguments of <5's in eq. (B.6). 
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